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ABSTRACT 



Non-thermal hard X-ray and high-energy (HE; 1 MeV < E < 100 GeV) y-ray emission in the direction of n Carinae has been recently 
detected using the INTEGRAL, AGILE and Fermi satellites. This emission has been interpreted either in the framework of particle 
acceleration in the colliding wind region between the two massive stars or in the very fast moving blast wave which originates in 
the historical 1843 "Great Eruption". Archival Chandra data has been reanalysed to search for signatures of particle acceleration in 
j) Carinae's blast wave. No shell-like structure could be detected in hard X-rays and a limit has been placed on the non-thermal X-ray 
emission from the shell. The time dependence of the target radiation field of the Homunculus is used to develop a single zone model 
for the blast wave. Attempting to reconcile the X-ray limit with the HE y-ray emission using this model leads to a very hard electron 
injection spectrum dN/dE oc E~ r with T < 1.8, harder than the canonical value expected from diffusive shock acceleration. 
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1. Introduction 

77 Carinae is a binary system composed of a massive pri- 
mary (M > 90 M Q , t] Car A) and a less ma ssive secondary 
(M < 30 M Q , 77 Car B) star (see e.g. iNielsen et all I20071) 
which orbit each other in 2022.7 + 1.3 days dDamineli et al.l 
• l2008al) . 77 Carinae experienced a historical outburst (the "Great 
£>T) ' Eruption") in the 19 th century and ejected ~12M o f gas which 
£f) \ move s outwards at an average speed of -650 km s _1 (Smit h~et~aTI 
I/") ■ 120031) forming the "Homunculus Nebula". Recent observations 
C*") ' show that 77 Carinae is surrounded at a distance of ~0.25pc 
* , by a very fast moving blast wave (3500 - 6000km s~') pro- 

■ duced in th e giant outburst of 1843 (also known as the "Great 
{SJ j Eruption") (Smith 2008). This blast wave currently overruns 
y— I . the "Outer Ejecta" a ring-like structure of material which orig- 

ILJ ' inates fro m an ejection of mas s from 77 Carinae ~ 500 - 1000 
. j-h ! years ago (Walbor n et al.ll978b . This fast-moving material mim - 

■ ics a low-energy supernova remnant (SNR) shell dSmithl 12008). 
$— 1 ' with a blast wave moving into the ISM wit h velocities com para- 
F$ I ble to the historical supernovae RCW 86 (IVink et al.ll2006l) and 

SN 1006 (IVinkll2005l) . 

There is evidence for the presence of relativistic particles in 
77 Carinae. The existence of non-thermal X-ray emiss ion was re- 
cently reported by the INTEGRAL collaboration dLevder et al.l 
2008). In the high energy (HE; 1 MeV < E < 300 GeV) domain, 
a source spatially coincident with the 77 Carinae position was 
repor t ed by the AGILE and F ermi collaborations (Tavani et al.l 
2009; lAbdo et~ai1l2009l l2010h . Within the measured positional 
uncertainties of INTEGRAL and Fermi-LAT, particle accelera- 
tion via the diffusive shock acceleration (DSA) process is pos- 
sible in the colliding wind region (CWR) of 77 Carinae and/or 
in the expandin g blast wave of the Great Eruption of 1843 (see 
lOhm et aLll2010i for a detailed modelling). The absence of sig- 



nificant variability in the 50 ke V and MeV-Ge V regime is sur- 
prising in a colliding wind binary (CWB) picture, particularly 
during periastro n passage where a collapse of the CWR is ex- 
pected (see e.g. iParkin et al.|[2.009l and references therein) and 
hence no particle acceleration should occur. The blast-wave sce- 
nario provides a good explanation of the observed emission be- 
cause of the existence of an extended emission region, thus ex- 
plaining the lack of significant variability of the source. While 
Fermi and INTEGRAL do not provide sufficient angular resolu- 
tion to resolve the blast wave, X-ray observations with the cur- 
rent generation of instruments should provide resolved images of 
the region. However there is likely to be some confusion of the 
possible non-ther mal X-ray shell with thermal emission from the 
Outer Ejecta (see ISeward et al.ll200ll for more details). In this 
work, archival Chandra X-ray data have been analysed to search 
for signatures of accelerated particles in the blast-wave region. 



2. Chandra data 

An 88.2 ks Chandra ACIS-I observation (obs ID 6402) of the 
Trumplerl6 region, taken in August 2006 was reprocessed us- 
ing the most recent version of the calibration files. Data pro- 
cessing and reduction was performed using the Chandra CIAO 
(v4.3) software package and CALDB (v4.4.3). The data were 
unaffected by soft-proton flares and so the full observation time 
was available for analysis. Standard procedures were followed 
and images (count-maps) were created in two energy bands; 
0.3 keV - 5 keV and 5 keV - 10 keV and are shown in Figure Q] 
The extended emission attributed to the expanding ejecta around 
77 Carinae can be seen clearly in the low-energy image. The 
high-energy image however is dominated by the bright point-like 
emission from 77 Carinae itself. Several other low-significance 
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Fig. 1. Chandra image of the region around 77 Carinae. The top 
(bottom) panel shows the counts map in the energy range 0.3- 
5.0 keV (5.0-10.0keV). The black ellipse and circle show the 
flux extraction region for the ejecta. See text for details. 



sources are detected above 5 keV, see iLevder et alJ (1201 Oh for 
details. 

Slices in Right Ascension and Declination were made 
through the low and high energy images and through a simu- 
lated image of the Chandra- ACIS PSF (created at the same chip- 
location as 77 Carinae in this observation). The slice regions were 
7" wide and the projections are shown in Figure [2] It is clear to 
see that there is little emission above the PSF in the high energy 
map (green points) from the 77 Carinae region. 

The flux from the shell was estimated from the elliptical 
region shown in Figure Q] A circular region (also shown in 
Figure [TJ centred on the position of 77 Carinae with a radius of 
0.305' was excluded from the flux calculation. The background 
was estimated from a large elliptical region to the south west 
of the "on region". All point sources were removed from the 
background region before extracting the flux. The background- 
subtracted flux from the shell in the 5-10keV band was calcu- 
lated to be 4.3xl0 _13 ergcm _2 s _1 . No spectral analysis of the re- 
gion has been attempted and so this value represents the upper 
limit on the non-thermal flux from the ejecta. 

3. Fermi-LAT data 

The hypothesis that particle acceleration is occurring in the ex- 



panding ejecta surrounding 77 Carinae suggested by Oh m et aD 



(120101) was based on an analysis of the first 1 1 months of Fermi- 
LAT observations. More recent analyses of increased data sets 
are now availab le and are used here to constrain the model of 
IOhmetal.l(l2010h . 

An analysis in cluding 21 m o nths o f Fermi-LAT data has 
been presented by iFarnier et al.1 d201 ll) . This analysis reveals 
two components of the HE emission; a low-energy part (here- 
after L-component) which is best described by a power-law 
with spectral index F = 1.69 + 0.12 and exponential cut-off 
at E c = 1.8+ 0.5 GeV, and a second, high-energy compo- 
nent (hereafter H-component) extending to « 100 GeV, best de- 
scribed by a pure power law with index 1.85 + 0.2 5. No tempo- 
ral va riability in either component is reported in IFarnier et al.l 
d2oTTh . H owever, some var i ability i n the H-com p onent is re- 
ported by I Walter & Farmer! (1201 ll) . IFarnier et all (1201 ll) have 
proposed that the Fermi-H component results from the interac- 
tion of accelerated protons and nuclei. This interpretation is at- 
tractive in the sense that accelerated protons can have a higher 
maximum acceleration energy and suffer less from losses than 
electrons. We follow this approach in this discussion attributing 
the two components to electrons (L-component) and protons (H- 
component) in the blast wave. 

4. Origin of the HE y-ray emission 

The properties of the HE y-ray emission from the region in and 
around r\ Carinae are extremely challenging to interpret within 
the framework of any current model. In particular, for models in 
which the emission originates in the CWR, the lack of observed 
variability for the bulk of the emission is hard to reconcile with 
the dramatic chang es seen at other wav elengths during perias- 
tron passage (see Dami neli et al.ll20"08bl and references therein) 
and t he very short cooling time of relativistic particle s in the sys- 
tem (IFarnier et al.ll201 It iBednare k & Pabich 2011). The outer 
blast-wave scenario proposed in lOhm et all d2010l) provides a 
promising alternative in the sense that short timescale variability 
is not expected. However, the Chandra observations presented 
here place rather tight constraints on this scenario with impor- 
tant consequences for shock acceleration in systems of this type. 
Here we present these constraints and the refined model of the 
blast wave emission and discuss the more general implications 
of our results. 



4. 1 . Improved Blast Wave Model 

The single-zone , time-dependent numeri cal model used in 
lOhm et all d2010l) and described in detail in lHinton & Aharonianl 
d2007l) has been modified here such that the radiation field en- 
ergy density of the Homunculus nebula, which had been as- 
sumed to be static is now modelled as a function of time. The 
Homunculus is illuminated from the inside by the star 77 Carinae 
with the (bolometric) luminosity L, r This starlight is reprocessed 
by dust in the optically thick nebula into IR light. In the steady- 
state case, its own luminosity Lh = L n , but shifted to the IR band. 
Its thermal time scale, t±, is given by f t h ~ Lth/Ln, with be- 
ing the thermal energy of the Homunculus nebula. For a mass 
of the Homunculus nebul a of 12 M and a temperature of 260 K 
E th dGehrz & Smithll 19991) is found to be about 10 45 erg. Using 
Lh ~ lO^ergs" 1 dCox etal.l ll995) leads to a thermal timescale 
of the order of a few weeks. This value is much smaller than the 
age of the Homunculus nebula and therefore, assuming a steady 
flux from the star 77 Carinae, a (pseudo-) steady state is a reason- 
able assumption. Any luminosity change of the star 77 Carinae 
would hence be followed by a corresponding luminosity change 
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Fig. 2. Slices through the Chandra images in Right Ascension 
(top panel) and Declination (bottom panel) (see text for details). 
The width of the extraction region was 7". Black points repre- 
sent the profile of the full energy range (0.3-10.0keV) image. 
Red points represent the low (<5keV) energy data and green 
points represent the high (>5 keV) energy data. The blue curve 
shows the simulated Chandra- AGS PSF at the chip-position of 
77 Carinae in this observation. The dashed vertical lines show 
the size of the excluded region around Eta Carinae (circular re- 
gion in Figure 1). The y-scale has been truncated to highlight 
the behaviour away from the central peak. Note that the PSF is 
normalised to the maximum of the full-energy-range data. 



4.2. Application to the data 

Fig.[3]shows the y-ray spectral energy distribution which is pro- 
duced by electrons accelerated in the blast wave and interacting 
with the time-dependent radiation field of the Homunculus neb- 
ula. The X-ray upper limits presented in this work restrict the 
model considerably. The spectral index of the accelerated elec- 
trons is constrained by the Chandra limit to be rather hard (r < 
1 .8). This value is much harder than the canonical value (F = 2) 
expected from diffusive shock acceleration, but could be realised 
in very strong shocks or in shocks which are modified e.g. by the 
pressure of the accelerated particl e population. Us ing a magnetic 
field strength of 10/uG (as used in lOhm etalll2010h . the total en- 
ergy in electrons E e for this model would be 6 x 10 45 erg, rep- 
resenting only a very small fraction (» 10~ 4 ) of the total kinetic 
energy Ey in the blast wave of w (4 - 10) x 10 49 erg which is in 
principle available for particle acceleration. For a slightly higher 
magnetic field strength of 20 /iG, and a spectral index of F = 2.1, 
the energy in electrons has to be smaller than 3 x 10 45 in order 
to agree with the X-ray limit. However, this would imply that 
the vast majority of the HE y-ray emission does not originate 
in the blast wave. Given these findings, the fraction of energy 
in non-thermal electrons compared to the total kinetic energy in 
the blast wave is even lower compared to the model presented 
before. The association of the soft y-r ay emission detect ed by 
INTEGRAL from the 77 Carinae region dLevder et al.ll2008l) with 
the blast wave is problematic due to the sharp, low energy cut- 
off required for c onsistency w i th the Chandra limits presented 
both here and in iLevder et al.l d2010l) . It seems likely that this 
ha rd X-ray emi s sion i s associated with the CWB as suggested 
bv ILevder et"ail ( l201 Oh rather than the blast wave. 



lOhm et aTl (120101) concluded that hadrons were less likely 
responsible for the single-component HE y-ray emission re- 
vealed in the first 1 1 months of Fermi-LAT data. This conclu- 
sion was derived from two factors; the presence of HE emission 
at « 200 MeV (below the =s 300 Me V threshold energy for n° 
production) and the fact that the maximum energy of protons as 
indicated by the curvature in the Fermi spectrum lay well be- 
low the expected maximum acceleration energy associated with 
either the age or size of the system. It has a lso been argued 
(lFarnieretal.1120111: iBednarek & Pabichll20TTh that the density 
of target material for pp interaction and subsequent 7r°-decay y- 
ray production would be too low in the 77 Carinae region. 



of the Homunculus. The temperature of the radiation field T 
at the location of the (expanding) ejecta is given by the Stefan 
Boltzmann law: 

T = (L H /(47r ( r4)) 1 

where, ^h(0 is the time-dependent radius of the Homunculus. 
Due to the mass of the expelled material Rn(i) follows free ex- 
pansion Rn(t) ~ v'h t , with v'h being the velocity of the ejecta of 
the Homunculus. Hence, the time dependent radiation field only 
depends on the evolution of the position of the blast wave with 
respect to 77 Carinae, the luminosity of 77 Carinae and the tem- 
pe rature of the radia t ion fie ld. Based on the historical light curve 
of iHumphrevs et ail (Il999h . and the fact that the IR emission is 
simply reprocessed star-light, we assume that the integrated lu- 
minosity of the nebula in IR is linearly rising with Eta Car visible 
magnitude my. 



The new Fermi-LAT data reveal the two-component na- 
ture of the HE emission, and present a good case that the H - 
component may be of hadronic origin dWalter&Farnierll2011h . 
The variability detected in the H-component point to an origin 
of at least part of this H-component flux in the colliding wind 
region. However it is likely that some part of this H-component 
originates from protons accelerated in the blast wave surround- 
ing 77 Carinae. Such emission would be non-variable and would 
contribute to the total flux of the H-component. Scaling the total 
flux in the H-component down by a factor of three (representing 
the decrease in flux in the high-energy Fermi component found 
bv lWarter&Farnierll201lL see green data points in Fig. [3} and 
using the available target material density of lOOcirT 3 leads to 
the requirement of 6 x 10 48 erg of energy in hadronic cosmic rays 
in the region. This represents 6 - 15% of the kinetic energy of 
the blast wave, given the uncertainties in the kinetic energy es ti- 
mates in the blast wave of 4-10 x 10 49 erg (ISmith et al.ll2003l) - 
conditions that can reasonably be met. 
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Fig. 3. Spectral energy distribution for the region within a few arc minutes of rj Carinae as described in lOhm et al.l (1201 Oh . Curves 
show a single zone time dependent model for continuous injection of electrons and protons (synchrotron in blue, inverse-Compton 
in red and 7r°-decay in green). A magnetic field strength of B — 10/vG and electron energy of E e = 6 x 10 45 erg is assumed for the 
model represented by the blue and red solid curves. A spectral index of F = 1.8 and maximum electron energy of £ malje = 1 10 GeV 
is used for this model. E p — 6 x 10 48 erg of energy in protons and a maximum energy of E maXtP = 2 TeV is required to reproduce 
the 7r°-decay y-ray co mponent indicated b y the green solid curve. Note that the green data points are simply the Fermi-LAT data 
points as presented bv lFarnier et al.l (1201 lh . scaled down by a factor of three. The dashed set of curves use a model which assume a 
magnetic field of B - 20 //G, E e = 3 x 10 45 erg of energy in electrons and spectral index r = 2.1 with the same maximum electron 
energy of £max,e = 1 10 GeV as used before. 



5. Summary and conclusions 

We have re-analysed archival Chandra data and placed limits 
on the non-thermal X-ray emission from the expanding ejecta 
surrounding ri Ca rinae. The single-zone numerical model of 
lOhm et all (1201 Oh has been adapted to account for the time- 
varying radiation field and to fit to recent Fermi-LAT HE y- 
ray data. The two-component nature of the HE emission is 
best explained by electrons (L-component) and protons (H- 
component) respectively. An attempt to reconcile the new limit 
on the non-thermal X-ray emission from the shell with the Fermi 
L-component data leads to a rather hard electron injection index 
of F <1.8. A steady hadronic emission component originating 
in the blast wave, at a flux level approximately one third of the 
total H-component, can be explained using reasonable numbers 
for the energy in cosmic rays, and the target density in the blast 
wave region and the maximum proton energy. 
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